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O Abstract 



Solar supergranulation remains a mystery in spite of decades of intensive studies. Most of the papers about supergranulation 
deal with its surface properties. Local helioseismology provides an opportunity to look below the surface and see the vertical 
structure of this convective structure. We present a concept of a (3+l)-D segmentation algorithm capable of recognising individual 
^ supergranules in a sequence of helioseismic 3-D flow maps. As an example, we applied this method to the state-of-the-art data 
and derived descriptive statistical properties of segmented supergranules - typical size of 20-30 Mm, characteristic lifetime of 
CO .18.7 hours, and estimated depth of 15-20 Mm. We present preliminary results obtained on the topic of the three-dimensional 
^structure and evolution of supergranulation. The method has a great potential in analysing the better data expected from the 
helioseismic inversions, which are being developed. 

Keywords: The Sun: photosphere. Sun: interior. Sun: helioseismology 
^ PACS: 96.60.Ly, 96.60.Mz, 96.60.Jw 
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1. Introduction 



Motions of plasma in the solar phot osphere have a few char - 



acteristic spatial scales (see review by iNordlund et al. ■ l2009h 
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Jlhe convective nature of the solar subsurface layers can be eas- 
ily seen in solar granulation, convective cells with character- 
istic scale of 1 Mm and lifetime of ~10 minutes. Exploding 
granules form a larger apparent structure, the mesogranulation. 
,Weak background magnetic field is not distributed randomly in 
the photosphere, but concentrated in lane-like structure on the 
boundaries of supergranules. Larger scales in motions extend 
further towards larger structures (giant cells, Rossby waves) 
iand end at a differential rotation. Understanding these motions 
is an important step towards understanding the solar dynamo, 
,a mechanism for origin and redistribution of solar magnetic 
fields. 

' It seems that the solar supergranulation is especially impor- 
tant in forming surface magnetic features. It has been known 



since 1950s from surface Doppler measurements ( e.g. lHartllI954 
Leighton et al.[ Il962t ISimon and Leightoni Il964l) . Since then. 



a vast amount of papers has been published with the topic of 
supergranulation. It is generally agreed that supergranules are 
cell-like structures with typical spatial scale of 30 Mm, life- 
time of ~24 hours, having predominantly horizontal velocity 
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field with RMS of ~300 ms and with basically no temper- 
ature contrast between the middle of the cell and its bound- 
ary. The origin of supergranulation is currently unknown. The 
classical mechanism invoked to explain the origin of supergran- 
ules is the latent heat of the recombinati on of He^^ into H e^ at 
roughly 10 Mm below the photosphere ( lGieraschlll985b . In a 
fluid at rest, such a heat release may trigger an instability, which 
can turn into motions at scales comparable to the supergranular 
scale. Recently, convection structures on supergranular scales 
have been se en in hydrodynamic s i mulations of subsurface con- 
vectio n (e.g. lRieutord et al.Ll2002tlStein et al. [ 120061: 1 lUstvugovl 



2008) 



The development of local helioseismology has made it pos- 
sible to investigate subsurface layers of the solar convection 
zone. Time-distance helioseismology relies on analysis of prop- 
agation of acoustic waves through the solar interior. The time, 
in which the wave packets travel from one point to another one, 
is affected by inhomogeneities along the propagation path. Us- 
ing various procedures, it is possible to invert the measured 
travel-time deviations in order to obtain information about the 
inhomogeneities in the solar interior, e.g. plasma flows. 

A traditional approach is an extension of the principles of 
geometrical optics, in which rays approximate wave-packets 
propagation. This approach has been used in many studies and 
it has been shown that it provides credible results for large-scale 
perturbations with a reasonable computational efficiency. New 
inversion methods using more realistic approximations (e.g. first 

order Born approximation) ar e being developed. Although the 

tests u sing simulated data (e.g. lZhao and Kosovichevll2003tlZhao et al 
2007h demonstrate the credibility of ray-based inversions for 
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the horizontal flow components in the first 5 Mm (and perhaps 
10 Mm), independent stud i es wi th a different inversion tech- 



cell. It has been shown in earlier works (e.g . lBeck and DuvaU Jr 
200ll;lGizon et al.Ll2003llHirzbergeret allboOSl) that a conve- 



nique (e.g. lJackiewicz et al.Ll2008l) suggest that already at 4 Mm 
depth the measured flow field is dominated by noise. This is- 



sue is currently being intensively worked on (e.g. lSvanda et al. 
in preparationi) . 



nient way to identify the individual convective cells, is the use 
of the horizontal divergence V/, ■ v. 



Flow maps computed by any method of local helioseismol- 
ogy record the structure and amplitude of flows in the certain 
layer of solar convection zone, given by the averaging kernel 
also computed by the inversion. The averaging kernel pro- 
vides important information about the smoothing of the result- 
ing flows. Another important parameter is an accuracy of the 
results, which is not precisely known in RLS ray-kernel based 
inversions (it is, howeve r, estimated as a re lative error of 10% 
in near-surface layers - iBirch et alj 1200 ll) . The accuracy is 



V . V = — -H — 

dx dy 



(1) 



Jackiewicz et al 



known in e.g. OLA-based inver sions (e.g 

2008 ; Svanda et al. , in preparationh . which take into account the 
covariance matrix of travel times and compute the propagation 
of errors in travel times into the inverted quantity. Helioseismic 
flow inversions allow for measurements of the flow maps at dif- 
ferent depths and therefore provide a 3-D tomographic image 
of the flows in the solar interior Therefore, they are suitable for 
studying the evolution of 3-D structure of convection structures, 
such as supergranulation. 

Three-dimensional inversions of travel times in quiet Sun 



where v i and v,. are horizontal components of the flow field. On 
the surface, V/, ■ v mimics the vertical component v^: the re- 
gions of a positive horizontal divergence are usually regions of 
up-flows, while regions with a negative horizontal divergence 
depict down-flows. It would be wiser to use v. directly in iden- 
tifying the supergranules, however it is not possible in current 
inversions, because a weak signal of v- is buried in noise and 
systematic errors, such as a crosstalk. 

A supergranule in a flow map near the solar surface is de- 
fined as a compact region of the positive divergence surrounded 
by the lanes of the negative divergence. In order to avoid the 
noise, thresholding and limitations on the detected structure 
size are necessary. The segmentation algorithm is, however, in- 
dependent of the supergranule definition and in principle should 
be usable for any reasonable definition. 
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regions have been shown in past studies, (e.g.puvall et al.L. , „ 
Kosovichev and DuvallL 1997 : Zhao and Kosovichev , 2003t" Duvall c 
20091) . The usual way of studying the vertical structure of su- 



pergranulation is the correlation bet ween the flow s at the given 
depth with the flows on the surface. iDuvalll ( 119981) showed that 
in the upper 5 Mm, the flows are correlated positively, then 
the correlation drops to zero and is negative in the range of 5- 
8 Mm, suggesting the existence of the "return flow". The cor- 
relation disappears below 8 Mm, which leads t o the depth of the 



superg ranulation of 8 Mm. Similar analysis bv lZhao and Kosoviche\ 



(|2003|) ended wi th the estimate of th e supergranulation depth 
of about 15 Mm. Svanda et al. I (l2009h compared the large-scale 
flow from helioseismic data-cubes with the surface flow derived 
from the tracking of supergranules in MDI Dopplergrams with 
a conclusion that the large-scale velocity field does not change 
much in upper 10 Mm depth and supergranules can be therefore 
treated as objects carried by the underlying flow. 

The main goal of our study is the development of an effi- 
cient code to study evolution of the three-dimensional structure 
of supergranulation. The paper is organised in a following way: 
in Section 2 we explain the segmentation method used to isolate 
individual 3-D supergranule and their histories in the sequence 
of helioseismic data-cubes. In Section 3 we briefly describe the 
helioseismic data we have at our disposal and the results we ob- 
tained after application of the described method. We conclude 
in Section 4 and discuss a future perspectives of the presented 
study. 

2. Segmentation method 

In order to isolate the supergranules and their histories in the 
sequence of 3-D flow maps we must first define a supergranular 
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Figure 1: An example of the horizontal divergence (positive part only, nor- 
malised to the maximum value) of the supergranular flow measured by time- 
distance helioseismology at four different depths. The data were smoothed by 
the low-pass Butterworth filter in order to suppress noise. We see that in the 
first three depths the pattern remains largely unchanged, i.e. depicting com- 
pact supergranular cells, while in the last displayed layer the pattern is dift'erent 
(slightly negatively correlated to the previous ones). The layer at the depth 
of 7.6 Mm is a transition layer between these two. The aim of the segmenta- 
tion routine is to find and connect the coiresponding areas with high value of 
positive horizontal divergence into compact three-dimensional object and track 
them in the consequent 3-D horizontal divergence data-cubes. 

In our case, we smooth all divergence maps in A:-space by 
a fourth-order low-pass Butterworth filter in order to supress 
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the structures that are significantly smaller than the expected 
size of supergranules. The threshold for the regions, which 
are considered as candidate supergranules, is selected as a e- 
multiple of the RMS-divergence at the given time and depth- 
layer As a result, the input of the segmentation routine is a 
set of consequent three-dimensional divergence maps, in which 
we keep smoothed regions with the horizontal divergence above 
the threshold, surrounded by the lane-like regions having zero 
value (cf. Fig. [1). The threshold parameter e is subject to the 
tuning, but remains constant for all times and depths, so that it 
takes into account possible systematic change in the flow am- 
plitude with time and depth. Here e- 1.5. 

The segmentation algorithm is based on the feature tracking 
approach, where we track supergranules from the surface down 
in depth and forward in time. In "(3+l)-D segmentantion" 
number 3 stands for the spatial dimensions and 1 for time. We 
start from the near-surface layer at the starting time fo, which 
corresponds to the beginning of the sequence of data-cubes. In 
the processed divergence map we label each continuous region 
with divergence values above the threshold by a unique tag - 
these are possible supergranules. These supergranular candi- 
dates are then tracked down to the deeper layers at the same 
time, assuming that the location of the same cell at dififerent 
depths does not change significantly. In this way, we select 
three-dimensional compact objects (3-D supergranules at time 
fo)- 

In the next step, we move to the consequent 3-D divergence 
data-cube at time t\ - + Af. We assume that during the 
time step Af, which divides two contiquous data-cubes (here 
Af = 8 hours), the supergranules did not move significantly, i.e. 
that supergranular candidates at a similar position represent the 
same supergranule at two consequent times. Given the difer- 
ential rotation rate in the solar photosphere, this assumption is 
very reasonable, the estimated shift is 2-A pixels with resolution 
of 1 .4 Mm. In the near-surface layer, we search for supergran- 
ules known from the previous time-step. If some of them are not 
identified, it is assumed that they are diminished and their his- 
tory ends here. Newly detected supergranules, without an ap- 
pearance in the previous time-step, are labeled with the unique 
tag. The search is then extended into the depth domain. As a re- 
sult, we obtain the first set of three-dimensional supergranules 
with parts of their histories, and a new set of three-dimensional 
supergranules, which do not have a history yet, because they 
have been seen for the first time. We then move to the conse- 
quent data-cubes at times f2 = fi -H Af, - f2 H- Af, . . . , and 
repeat the algorithm again, until we reach the end of the contin- 
uous sequence of data-cubes. The procedure results in a set of 
three-dimensional supergranules and their histories, describing 
their location in space and time. 

Physically it is possible that the supergranule splits in two, 
or that two supergranules merge into one. If this happens, the 
split or merged supergranules are labeled as completely new 
ones, ending the life of their predecesors. Supergranules not 
having a complete history, i.e. being on the edge of the data- 
cube, are removed from the sample. 
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Figure 2: A comparison of the real measured Dopplergram and the synthetic 
one, constructed from the time-distance data. We see that the match is not 
perfect, however, still reasonable. Correlation coefficient in this example is 
0.37. 



3. Results 

3.1. Data 

We have developed and tested the method on data-sets avail- 
able today, which consist of a series of 3-D flow data-cubes 
from local helioseismology (see example in Fig.[T]). These flow 
data-cubes were obtained by r ay-kernel based time-distance 



helios eismology dDuvall et al.L II993I: iKosovichev and Duvall , 



19971) using phase-speed filters and RLS fitting, and are con- 
sidered state-of-the-art. The data we have at our disposal w ere 
used in previous studies (e.g. IKosovichev and Duval 



Svanda et al. 



2009h and consist of 46 flow-map sets in the vicinity of the ac- 
tive region NOAA 9393/9433. This sequence of three-dimensional 
flow maps was recorded in March and April 2001. To obtain 
reasonable signal-to-noise ratio, each of these maps is averaged 
over 8 hours and also in spatial coordinates, so that the eff'ective 
resolution of the results is 5.6 Mm. The data-cube contain two 
temporal discontinuities, when the active region was situated on 
the far side of the Sun. 

The data were originally computed for studies of the depth 
structure and evolution of a large active region. The presence of 
the magnetic field aff'ects convection and also the propagation of 
seismic waves, therefore magnetic regions are not suitable for 
studies of supergranulation. For our purpose, we use only quiet 
Sun regions around the magnetic areas. This selection limits the 
effective field-of-view and also the number of observed super- 
granules that can be used for statistical studies. As we intend to 
develop the method of (3+l)-D segmentation on available data 
and use it for serious research using better data available in the 
future, we do not see the relatively small amount of supergran- 
ules as a serious drawback. 

3.1.1. Data testing 

In order to test the credibiUty of the time-distance data we 
have at our disposal, we constructed the synthetic surface Dopp- 
lergrams and compared them with the observed ones. One has 
to keep in mind that we do not compare exactly the same fea- 
tures: synthetic Dopplergrams are constructed by taking the 
line-of-sight projection of the horizontal (i.e. Vv/v,,) velocity 
field representing flows at the depth of 0.77-1.77 Mm. We do 
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Figure 3: A histogram of lifetimes of supergranules detected by tlie segmen- 
tation method. A fit following an exponential decay is over-plotted with an 
expectation lifetime of 18.7 hours. 



not have pure surface helioseismic measurements at our dis- 
posal, because the surface /-mode is filtered out during the 
application of time-distance helioseismology used for the in- 
ve rsion. Moreove r, vertical v. component is unreliable to use 



jZhao et al .1120071) . Therefore, we cannot expect a perfect match 
between the reconstructed and measured Dopplergrams. 

Nevertheless, as we see in Fig.|2] we find many similarities. 
The correlation coefficient is typically 0.4 and it is higher at 
around 55 degrees from the disc centre. Based on the direct 
comparison with the surface measurements, we may estimate 
that the accuracy of helioseismic data lies between 10 and 20% 
of relative error 

3.2. (3+1 )-D segmentation 

As already mentioned, we define the supergranule as a com- 
pact region of the positive horizontal divergence in the near- 
surface layers. Therefore, we first use horizontal components of 
the measured plasma flow for computation of the horizontal di- 
vergence. Resulting (3 +l)-dimensional divergence data-cubes 
are smoothed by a fourth-order low-pass Butterworth filter with 
a cut-off of 0.034 Mm ' in order to suppress the random noise 
in the data and enhance the structures on supergranular scales. 

Following that, the segmentation algorithm described in Sec- 
tion 2 is applied. We recognised 1212 supergranular cells and 
their histories. In the sample, short living cells having volume 
less then 300 Mm-^ are designated as noise and removed from 
the set. This value is roughly one tenth of the expected volume 
of supergranular cell having 30 Mm in diameter and depth of 
5 Mm. 

We investigated the resulting set of supergranules to confirm 
that its statistical properties are consistent with hterature. 

A histogram of the lifetimes of supergranules in the sample 
is displayed in Fig.|3] An exponential fit of the form 



f(T) ~ exp 



(2) 



(120041) . From our analysis, it becomes clear that supergranules 
with higher values of the horizontal divergenc e in near-surface 
layers depict longer lifetime, in agreement with lHirzberger et al 
(120081) . The dependence is almost linear for supergranules hav- 
ing maximum horizontal divergence values 10-60 jU5 for su- 
pergranules having maximum horizontal divergence larger than 
60 jU5"' the lifetime rem ains almost constant (~30 hours), again 
with an agreement with iHirzberger et al.l ( 120081) . The typical 
size of the detected supergranular cells is 20-30 Mm. The life- 
time of the cells with this size is around 25 hours. 

Having supergranules defined as three-dimensional objects 
in the time sequence makes it possible to study their evolution. 
We do not observe significant rise or decline of the typical su- 
pergranule during its lifetime. The gravity centres of the hori- 
zontal divergence regions remain at the depth of 3.2 Mm over 
almost the whole life-span. We measure an insignificant aver- 
age rise of the gravity centre in the growing phase of the su- 
pergranules by roughly 0.2 Mm with the speed of 3 m s"' , both 
values are beyond the precision of the used data. In some super- 
granular histories, we detected a initial rise of the region with 
the maximum divergence with the speed of (300+100) ms"'. 
These observations unfortunately strongly depend on the qual- 
ity of the used data, which might be in doubt in our case. There- 
fore in the present analysis, we are far from making conclusion 
about the physics of supergranules. 

The evolution of the average supergranule in time-depth 
space is displayed in Fig. |4] This is a result of averaging over 
215 supergranules which survived for more than 24 hours. All 
supergranules were co-aligned in time to the moment, when the 
maximum of the divergence in the near-surface layer of 0.77- 
1.77 Mm were established. The characteristic horizontal diver- 
gence used in this plot was calculated as a horizontal average of 
the horizontal divergence in the supergranular cell at the given 
time and depth. 

The overall picture on the plot again depends highly on the 
accuracy of the helioseismic data, especially whether the data 
are free of systematic errors. Nevertheless, current results con- 
firm that supergranules are of convective origin. They consist 
of the region with positive divergence in the sub-surface layers 
and the region of the negative divergence, the region of the "re- 
turn flow", extending deeper down. The upper "outflow" region 
evolves slowly in the intensity over the course of supergranular 
life, but remains confined in the depth till ~6 Mm. The "inflow" 
region seems to extend almost twice as deep at the maximum 
phase than during the growing and decaying phases. Based on 
the displayed results we estimate the depth of supergranulation 
of ~15 Mm at the growing and decaying phases, and ~20 Mm 
at the maximum phase. 

With a certain amount of phantasy we could say that the 
average supergranule starts to form 50 hours before its maxi- 
mum at the depth of ~7 Mm and extends towards the surface 
and deeper down as the time progresses. The evolution in time 
is slightly asymmetrical showing insignificantly slower decay 
compared to the formation. In the average picture, we do not 
see signs of structural rise or decline in depth. 



provides an expect ation lifetime of t = 18.7 ho ur s, which is in 



agreement with e.g. lDeRosa and Toomrd (120041) or lDel Moro et al. 
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Figure 4: The average supergranule evolution in time and deptli, averaged over 
215 supergranules with the lifetime longer than 24 hours. Shades of grey depict 
the average value of the horizontal divergence in the supergranular cell in the 
given time and depth. 



4. Conclusions 

We have developed and tested the method of (3+l)-D seg- 
mentation of supergranules as evolving three-dimensional ob- 
jects in the series of helioseismic three-dimensional flow maps. 
The segmented supergranules display surface statistical proper- 
ties that are in agreement with widely accepted parameters, e.g., 
typical size of 20-30 Mm, characteristic lifetime of 18.7 hours, 
and estimated depth of 15-20 Mm. We made an attempt to 
investigate the evolution of the average supergranule in time 
and depth, which could confirm the convective nature of super- 
granules. The results do not a ddress the issue of tendency of 
supergranules to support waves dGizon et aUl2003h . The wave- 
like properties manifest as apparent motions of supergranular 
structure, which are not decomposed when tracking individual 
supergranular cells in space during the process of segmentation. 

From our analysis it seems that supergranules extend deeper 
at the maximum phase of their evolution than during the phases 
of growth and decay and have slightly asymmetric evolution in 
time with decay slower than growth. 

These results are, unfortunately, sensitive to the quality of 
the data used for segmentation. We are well aware of papers 
which doubt ray-kernel based inversions showing that wave- 
effects have to be taken into account (e.g. iJensen et al. ■ I2OOOI: 
Birch and Kosovichev . l2000h as the wavelength of seismic waves 
is comparable with the size of inhomogeneities in the Sun. It 
is therefore likely that wave-based kernels could yield depth 
inversion results that are qualitatively different from the ones 
based on ray kernels. However, the data used in this study are 
still considered state-of-the-art and provide us with sequences 
of data-cubes suitable for developing and testing of the described 
method. Born-kernel based time-distance inversions for flows 



are being deve loped (e.g. lJackiewicz et al.Ll2008tlSvanda et al 



in preparatiorj) . especially in the frame of observations com- 
ing from Solar Dynamics Observatory. It will still take some 
time before the higher-quality flow maps suitable for (3h-1)-D 
segmentation will be available. We are looking forward to this 



data and believe that our method applied to these data-sets will 
provide a deeper insight into the 3-D structure and evolution of 
supergranulation. 
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